Persistence of Li Induced Kondo Moments in the Superconducting State of Cuprates 
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Using 7 Li NMR shift data, the anomalous local moment induced by spinless Li impurities persists 
below T c in YBa2Cu30e+j/. In the underdoped regime, the moments retain their Curie law below 
T c . In contrast, near optimal doping, the large Kondo screening observed above T c (Tk = 135 K) 
is strongly reduced below T c as expected theoretically when the superconducting gap develops. The 
limited spatial extent of the induced moment (on 1 st near neighbour Cu) is not drastically modified 
below T c , which allows a comparison with STM determination of the local density of states. Our 
results constrain theoretical models of the impurity electronic properties. 



The influence of impurities on superconductors has al- 
ways been used as an effective probe of their actual prop- 
erties. For example, while magnetic impurities are the 
most prominent s-wave pair breakers, any type of scat- 
tering is detrimental to d-wave superconductivity. To our 
knowledge, experimental investigation of the modifica- 
tions of the magnetic properties of an impurity below the 
superconducting transition T c has never been performed. 
Macroscopic bulk magnetic experiments are unadapted 
since the various contributions to the susceptibility can- 
not easily be singled out below T c . Local measurements of 
the susceptibility of these moments using hyperfine tech- 
niques are in principle possible, but are usually prohib- 
ited by technical limitations, such as strong spin lattice 
relaxation effects for the impurity NMR. 

Despite this experimental void, an extensive theoret- 
ical work has been devoted to this question in classical 
Fermi liquids Q Q . The behavior of the moment below 
T c is predicted to depend both on the shape of the su- 
perconducting gap and on the Kondo temperature Tk of 
the moment in the normal state. Tk is a signature of the 
screening of the moment by the conduction electrons and 
is related to its coupling J to the carriers. The primary 
effect, anticipated, but not observed so far, is a reduction 
of the Kondo screening due to the pairing of the carriers. 
For small J, this results in a complete restoration at low 
T of the Curie susceptibility of the moment. 

In cuprate superconductors, which are correlated elec- 
tronic systems, the magnetic properties of impurities are 
more intricate. In the normal state, spinless impurities 
like Zn g {§ or Li @ g substituted on the copper site 
of the Cu02 planes induce a local moment in their vicin- 
ity. This moment extends essentially on the four Cu near 
neighbours (n.n.) of the impurity. Its static J5| and dy- 
namic |(| susceptibilities exhibit a Kondo like behaviour 
with a large range of Tk values, which can be spanned 
by changing hole doping. The effect of superconductivity 
on this moment addresses two issues: the persistence of 
magnetic correlations in the superconducting state, and 



the influence of d-wave pairing on the Kondo screening. 

We present here the first measurements of the induced 
moment properties below T c . they are performed using 
7 Li NMR since the transferred hyperfine couplings are 
weak enough that relaxation effects do not prohibit NMR 
spectroscopy of the moment. We propose a method to 
extract the susceptibility xioc probed at Li sites. This re- 
quires correcting the internal field seen by Li for screen- 
ing and vortex effects due to superconductivity. It will 
then be shown that the Li induced moments survive be- 
low T c and that Tk is strongly reduced. Furthermore, 
we will demonstrate that these induced moments remain 
confined primarily to the I s * n.n. coppers below T c . Re- 
cent scanning tunneling microscopy (STM) experiments 
in the Zn substituted Bi2212 cuprate gave a measure of 
the local density of states (LDOS) in the superconduct- 
ing state [0. The STM data suggest the occurence of a 
LDOS peak near the Fermi level on the Zn site and on 
the 2 nd n.n. Cu. This location contrasts with our finding 
of a magnetic state located dominantly on the 1 st n.n. 
Cu. The discussion of this discrepancy will lead us to 
favor theoretical models which incorporate the magnetic 
character of the impurity. 

The Li substituted samples YBa2(Cui_ Xn Li Xji )306+j / 
are those used in (5|. The two batches with Li nominal 
concentrations x n = 1% and 2% had an effective in-plane 
Li concentration of 0.85% and 1.86% per CuC>2 layer. 
Two oxygen contents were obtained from each batch cor- 
responding to optimally {y = 0.97) and under- (y = 0.6) 
doped regimes. Their T c were found to be 85.3 K and 
79.5 K at optimal doping, and 41 K and 25 K for the un- 
derdoped materials. The sample crystallites were aligned 
along the c crystallographic axis with an applied field. 
This allows accurate NMR measurements performed in 
fields parallel to c ranging from 3 to 7 Tesla. Below 
T c , NMR spectra are too broad for Fourier Transform 
spectroscopy and were then measured point per point by 
sweeping the frequency over a few hundreds of kHz 0] . 

In the superconducting state, the internal field B at 
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any point in the sample may differ from the applied field 
B app due to large screening effects. NMR measurements 
of the spin contribution to the NMR shift, which is pro- 
portional to xioc are not straightforward. The spectral 
position w*of the NMR signal of a 7 Li nucleus is related 
to its NMR shift 7 K by 



j(B app + SB)(l+ 7 K) 



range of temperatures below T c JT3] 14 1. Each vortex is 
then moving much faster than the NMR timescale, which 
averages out both the broadening and the screening ef- 
fects, leading to SB = 0. The exact determination of the 
melting temperature, which depends on B appi well 
as on the sample microstructure, is beyond the scope of 
this work. In the underdoped compound, for B app = 7 
(1) T and T > 10 K, the vortices should always be in the 



with SB = B — B app , and 7 the nuclear gyromagnetic 
factor. An independent determination of the distribution 
of SB at Li sites is needed. We can use the fact that the 
average of SB is almost independent of B app in a large 
range of fields. Indeed, the magnetization M oc SB has 
been measured to be flat from B app = 2 to 12 Tesla in 
similar pure powder ceramics [ [Tcfl . Measurements for two 
values of B app yield both SB and K at the Li sites from 
Equjl] applied to the peak position of the NMR. 

Let us specify first how lo* has been extracted from the 
NMR spectra such as those plotted in Fig.l. They con- 
sist of three lines at lo* — lo c , lo* , and lo* + lo c , due to the 
quadrupolar splitting of the 7 = 3/2 Zeeman transitions 
by the electric field gradient (EFG) at the Li site. The 
quadrupolar frequency lo c is proportional to the EFG in 
the direction c of B app ||. Above T c , the two outer lines 
are broader than the central one due to a small distribu- 
tion of lo c (typically Slo c Ilo c ~ 10%). The other source of 
broadening, common to the three lines is the local distri- 
bution of hole content which induces a slight distribution 
of 7 K. It scales with 7 K and therefore increases with de- 
creasing T. Below T c the presence of pinned vortices 
induces a distribution of SB amoung Li nuclei, leading 
to an additionnal asymmetric broadening (Eq.^), as ob- 
served in Fig.l. The high frequency tail and the cen- 
ter peak of the line correspond respectively to Li sites 
in the vortex cores with SB > 0, and between vortices 
with SB < f§]. We fitted the T <T C spectra using the 
T > T c shape convoluted by an asymmetric gaussian rep- 
resenting the SB distribution. The resulting values of SB 
using two measurements of the central line peak position 
lo* are plotted in Fig. 2 for optimal doping. The nega- 
tive sign of SB confirms that the susceptibility probed by 
the peak position lo* is associated to Li defects between 
vortices in the bulk superconducting state. At low T, the 
obtained SB ~ —30 G is consistent with measurements 
by ^SR or 89 Y NMR in pure compounds Q. The T- 
dependence of SB originates from the T variation of the 
superconducting screening and of the field distribution in 
the vortex network. The data for 8B(T) has been fitted 
by a phenomenological power law corresponding respec- 
tively for x n = 1% and 2% to SB 1% = -32(1 - (T/77) 2 ) 
and SB 2 % = -26(1 - (T/60) 15 ) in Gauss units. The 
decrease of SB with Li content can be explained by the 
concomitent increase of the penetration depth A mea- 
sured by fiSR JlI) . Above 77 and 60K respectively, we 
found SB = within error bars (±4 G). For optimally 
doped samples and applied fields of a few Tesla, it has 
been shown that the vortices are in a liquid state in our 
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lq ]. Indeed we find that SB = within 
accuracy, so that no correction to lo* was 
needed for such conditions. 

Using this determination of SB in Eq.|l|, the shift 7 K 
can then be safely extracted. As seen in Fig.l, this 
shift and therefore xioc are increasing with decreasing 
T below T c . This is confirmed by systematic measure- 
ments of 7 K for all concentrations of Li and oxygen 
dopings represented in Fig. 3. A more compelling repre- 
sentation of the variations of 7 K is obtained by plot- 
ting 1/{ 7 K — 7 Kq) ~ l/xioc versus T as done in Fig. 4. 
7 Kq is the T independent part of the shift and is mea- 
sured from high T data to be much smaller than the 
observed variations of 7 K . In this plot the Curie- Weiss 
law 7 K — 7 K Q — C/(T + Tk), which represents xioc in 
the normal state is a straight line with slope C _1 which 
intercepts the horizontal axis at —Tk- 

In the underdoped regime, from both Fig. 3 and 4, no 
significant change occurs at T c , i.e the almost perfect 
Curie law observed above T — 80 K is not affected by 
superconductivity and corresponds to Tk = 2.8 ± 1 K 
|16| . To our knowledge, this is the first measurement of 
Xioc in a superconductor for moments appreciably cou- 
pled to the carriers. In contrast we find a sharp increase 
of 

Xioc below T c at optimal doping, as seen in the lower 
panel of Fig. 3. This increase is not expected for a usual 
Kondo impurity. This can be seen by comparing our 
data with a prototype of the Kondo susceptibility such 
as XFe of Fe impurities in the dilute alloy CuFe. We 
have scaled the T variation of XFe measured in fll7| by 
a factor T K (0 7 )/T K (CuFe) ~ 135/27.6 ~ 4.9 to fit the 
normal state data. One can obviously see in Fig. 3 and 
4 that the rescaled XFe saturates at low T like xioc but 
at a much lower value (of about a factor three). The 
data for xioc(T < T c ) is better fitted by another scal- 
ing of XFe also represented in Fig. 3 and 4, leading to 
Tk = 41 ± 7 K instead of 135 K above T c . The moments 
survive below T c even at optimal doping. They still dis- 
play a Kondo-like susceptibility with a weaker screening 
than in the normal state. This value of Tk is consistent 
with the analysis of the specific heat measurements of 
Zn substituted YBaCuOy by Sisson et al. |Q who at- 
tributed the absence of a Shottky anomaly below T c to 
Kondo screening of the Zn induced moments. 

In the superconducting state of a Fermi Liquid, the de- 
crease of the DOS at the Fermi level prohibits the devel- 
opment of the Kondo divergence near or below a critical 
coupling J c . This is also true for d-wave superconductors 
for which the gap corresponds to a linear energy depen- 
dence of the DOS Q j|, except at low T in the presence 
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of impurities. Qualitatively this explains the reduction 
of Kondo screening and Tk at optimal doping below T c . 
Renormalization group numerical studies using a realis- 
tic set of parameters may quantitatively account for the 
behavior observed in Fig. 4 at optimal doping 0. 

In the underdoped regime, the absence of detectable 
modification of xioc below T c could result from the al- 
ready small value of Tk found in the normal state. In 
contrast with optimal doping, any reduction of Tk be- 
low T c cannot be observed in our experimental condi- 
tions where T 3> Tk- The low value of Tk already in 
the normal state could be explained by the occurence 
of the pseudogap, similar to the effect of the supercon- 
ducting gap at optimal doping. Both gaps display the 
same d-wave symmetry |l9| , and should lead to a reduc- 
tion of Tk- However, the significant difference between 
the optimal and underdoped regimes in the apparent 
Tk (respectively 41 K and 2.8 K) in the superconducting 
regime remain to be understood. 

Even though the above Fermi liquid picture explains 
the T behavior of xioc it cannot account for the very 
existence of the moment induced by spinless impurities. 
This moment is the result of electronic correlations in- 
trinsic to the pure cuprate. From NMR experiments per- 
formed in the normal state, this moment consists of a 
staggered AF state which extends on many lattice sites, 
but resides predominantly on the impurity 1 st n.n Cu 
sites. The present experiment demonstrates that this is 
still true below T c . At optimal doping, scaling with \Fe 
yields numerical values for the Curie term C of 5.1 10 4 
K.ppm and 4.8 ± 1.1 10 4 K.ppm above and below T c re- 
spectively. These values are only 25% smaller than in the 
underdoped case. Hence C is almost unaffected either by 
superconductivity or hole doping. Therefore the interac- 
tion with the charge carriers does not modify the effective 
moment on the 1 st n.n. but merely induces a modification 
of Kondo screening. 
The measurements by STM on Zn substituted Bi2212 
also imply a small spatial extent of the impurity LDOS 
below T c . This LDOS exhibits a narrow resonance peak 
at an energy of 1.5 meV («18K) below the Fermi level. 
This energy scale is close to our Tk value, suggesting a 
common origin for the two phenomena. But the absence 
of LDOS on the 1 st n.n. found by STM contrasts with our 
finding. Most computations of the LDOS can reproduce 
the symmetry of the observed STM pattern, but did not 
consider the existence of an induced magnetism |2(| [^l) . 
In order to explain altogether the magnetism and the 
STM experiments, we suggest two possible scenarii: 

i) For an on site potential, some interpretations |2C| l 
predict a weak LDOS on the impurity site, and an en- 
hanced LDOS on the 1 st n.n., which might agree with 
NMR. If, as usually assumed, the tunneling is vertical 
from the Bi to the Zn site, this prediction would result 
in an observation of a large LDOS on the 1 st n.n., in 
contrast with STM data. We suggest that the tunneling 
occuring through a BiO layer could have a much larger 
matrix clement to the n.n. sites than to the one on the 



vertical. This would yield an apparent LDOS on-site and 
on the 2 nd n.n in the STM experiments H). 

ii) Let us assume that the perturbation induced by the 
impurity corresponds to a potential on its n.n. Cu. We 
can anticipate that calculations along the lines of p0[ , 
with a potential on a n.n. Cu, would induce an LDOS 
on the impurity and the 2 nd n.n. Cu sites. The total 
actual LDOS would then be consistent with STM for a 
vertical tunneling. Polkovnikov et al. have done such a 
computation |^2| , which is furthermore the first realistic 
attempt to account for the magnetic properties detected 
by NMR. They introduced an extended moment on the 
four Cu n.n. sites, exchange coupled to the quasiparticles 
of the superconductor and could reproduce the spatial 
dependence of the LDOS observed by STM. 

In conclusion our measurements show that the moment 
induced by spinless Li in the normal state of cuprates still 
exists in the superconducting state. When the Kondo 
temperature Tk is comparable toT c , the Kondo screen- 
ing is strongly reduced below T c . This is consistent with 
computations made in a classical Fermi liquid in the pres- 
ence of a d-wave gap. This reinforces the analogy of the 
magnetic behaviour with that of a classical Kondo effect. 

We have found that the dominant magnetic contribu- 
tion still resides on the 1 st n.n. Cu of the impurity below 
T c . Therefore the short range AF correlations remain in 
the superconducting state. The present work leads us 
to believe in a common understanding of the local mag- 
netism and LDOS, in the spirit of p2[ . Determination 
of the evolution of the LDOS with temperature and hole 
doping should help to establish the relationship between 
NMR and STM results. This would also constrain the 
theoretical models, which should in addition account for 
the low T magnetic susceptibility using the energy de- 
pendence of the LDOS. 
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FIG. 4. T variation of the inverse of 7 K(T) - K which 
represents the inverse of the local susceptibility xiocof the in- 
duced moments nearby Li. Solid triangles correspond to un- 
derdoped x n — 1% YBa2Cu306+ H while solid (empty) circles 
correspond to x n — 1% (2%) at optimal doping. The fits with 
the CuFe Kondo susceptibility of fig. 3 are plotted with the 
same symbols. 



FIG. 1. 7 Li NMR spectra for YBa 2 Cu 3 6 .97 with x n = 1% 
obtained either by Fourier Transform or point by point. In the 
superconducting state 7 K has been obtained after correction 
of demagnetization effects using EqEI. 



FIG. 2. Difference SB between applied and internal field 
on Li sites in optimally doped compounds deduced from mea- 
surements in B apP = 3 and 7 Tesla. Straight and dot lines are 
phenomenological fits given in the text. 



FIG. 3. T variation of the 7 Li NMR shift 7 K . The arrows 
indicate T c in a 30 G applied field. For the underdoped 06.6 
sample, the full line is a Curie-Weiss fit with C = 6.5 10 4 
ppm.K and Tk = 2.8 K. For the optimally doped O6.97 sam- 
ples , the prototype Kondo susceptibility of Fe in CuFe from 
Ref. jl7| scaled to fit our data above (below) T c is represented 
by a full (dashed) line. 
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Figure 1 (Bobroff et al.) 
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Figure 2 (Bobroff et al.) 
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Figure 3 (Bobroff et al.) 




Figure 4 (Bobroff et al.) 



